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1 Introduction 



In the modern age of nanoscience and technology, spin dependent electron transport has opened 
a new branch both in theoretical and experimental Condensed Matter Physics, the so-called 
Spintronics. The word spintronics (spin transport electronics) was originated in 1994 during the 
manifestation of a new electronic device based on the role of the spin of the carriers, instead 
of or in addition to their charge. The advanced technological progress in fabrication and char- 
acterization of hybrid nanostructures has enabled us to study the spin dependent transport in 
low-dimensional systems like quantum wires, quantum dots, array of quantum dots, etc., in a 
very tunable environment. 

For modeling and building spin based electronic devices, a clear understanding of spin depen- 
dent transport through hetero-interfaces is needed. Presently, spintronic devices are manufac- 
tured with the aid of two different approaches. The first approach is associated with the develop- 
ing of the existing giant-magnetoresistance (GMR)-based technology. A giant-magnetoresistive 
sandwich structure is a prototype device which has widespread application in industry as a 
read head and a memory-storage cell. This structure consists of alternating ferromagnetic and 
nonmagnetic metallic layers, and, changing the relative orientation of the magnetization vec- 
tors in different magnetic layers one can tune the device resistance from the very low value to 
a very large one. This phenomenon can be utilized to sense changes in magnetic fields. The 
second approach searches novel techniques for generation as well as utilization of spin-polarized 
currents. The investigation of spin transport in semiconducting heterostructures is a major chal- 
lenge to us, and, people are searching the ways in which these systems can act as spin polarizers 
and spin valves. The reason is that the existing metal-based devices cannot amplify signals, 
while semiconductor-based spintronic devices can in principle provide amplification and serve 
as multi-functional devices. Not only that, the semiconductor-based devices can be integrated 
much easily with the traditional semiconductor technology. 

At present the generation of spin polarized electron beam and its transmission through mat- 
ter is one of the most promising research topics in spintronics. At the same footing, several 
other phenomena have to be explored to develop the spin dependent transport. For our illus- 
trative purposes, here we mention some of them as follow. In metals, the excitation as well 
as the switching of magnetization with the help of spin torques is an important issue. The 
spin dependent transport in metal-semiconductor interface, the spin tunneling through insula- 
tors, semiconducting heterostructures, etc., also draw much attention in this particular field. 
All these studies may be helpful for the development of quantum information processing and 
quantum computation. 

In addition to the spin dependent transport measurements, the studies of spin polarized 
beam support an interesting model, the so-called 'two current model', which can be used for 
those electrons having kinetic energies ranging from diffusive to ballistic regime. The model 
illustrates that, electrons with up and down spins form two different current carrying channels. 
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associated with the different resistivities. This phenomenon provides dramatic effects, when the 
magnetization changes abruptly at the interface between two magnetically distinct materials. 
As representative examples: GMR Effect, TMR Effect (tunnel-magnetoresistance). Spin Torque 
effect during spin injection and spin transport across metal-semiconductor interfaces. Let us 
now concentrate our study on the generation of spin polarized electron beams. 

2 Generation of spin polarized electron beams 

2.1 Separation of electronic spin states 

To reveal the interaction of spin polarized electrons with matter, we need a source of spin 
polarized electrons [1]. In the pioneer experiment of Stern and Gerlach, the two spin states 
of the unpaired valence electron of the charge neutral Ag atom were separated by applying 
an inhomogeneous magnetic field. However, for a free electron there is no such method which 
can separate the two spin states by applying an external magnetic field. This is due to the 
overwhelming effect of the Lorentz force acting on the charge of moving electrons. 

It is well known that ferromagnets are the natural sources of polarized electrons, and one 
can extract the polarized electrons using photo emission technique or through proper tunneling 
method. But, there are some disadvantages for the ferromagnetic sources, since to switch over 
the spin direction, the inversion of the magnetization of the cathode is needed. This can make 
a trouble in emerging electron beam due to the existence of the residual stray magnetic fields 
at the ferromagnetic cathode. 

In the present age, the widely used source of polarized electron is the Ga-As type semiconduc- 
tors, where electrons are extracted through photo emission. Though Ga-As is a nonmagnetic 
semiconductor, but a large amount of spin polarized electrons can be extracted with the aid of 
a circularly polarized infra-red radiation. 

2.2 Ga-As semiconductor as a polarized electron source 

Here we explore how a Ga-As semiconductor can be used as a source of polarized electrons. 
In Ga-As material, the electrons at the top of the completely filled valance band are localized 
between the As atomic sites (Z = 33) and the electrons at the bottom of the conduction band 
are localized on the sites of Ga atoms {Z = 31). Accordingly, the electrons can undergo direct 
vertical transitions from the top of the valance band to the bottom of the conduction band. 

The actual mechanism of getting spin polarized electrons from Ga-As source is schematically 
shown in Fig. [TJ The valance band of the Ga-As semiconductor is made up of 4p atomic orbitals 
of As atoms those are splitted by spin orbit (SO) interaction into 4p|/2 ^^'i ^P3/2 states with an 
energy separation ~ 0.34 eV. 4pi/2 state is not shown here in the figure (Fig. [1]). On the other 
hand, the conduction band is composed of 4s atomic orbitals of Ga atoms giving rise to a 4s^/2 
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state. By choosing the photon energy ~ 1.5 eV (~ band gap energy), higher energy transition 
i.e., the transition from the state to the state 4si/2 can be blocked. 

Due to the existence of the cubic crystalhne symmetry in the electric field, the 4p3/2 state 
is four-fold degenerate in energy. The conservation of angular momentum in the absorption of 
circularly polarized light provides Arrij = ±1 for the allowed transitions. With the help of left 
circularly polarized (LCP) light, transitions with Arrij = — 1 are allowed only, while, for the 
transitions with Am^ = +1, the right circularly polarized (RCP) light is needed. Thus, in cubic 
Ga-As only the transitions (I) ^Pz/2{f^j = +3/2) 4si/2("^j = +1/2) and (II) 4p3/2(w.j = 1/2) 
— > Asii2{mj = —1/2) can occur with the help of LCP light, obeying the selection rule Am^ = — 1. 
A similar type of transitions also take place from the other two states in the valance band to 
the two states in the conduction band by using the RCP light, following the selection rule 
Am,- = +1. Therefore, both the spin states will be present in the conduction band. But, the 
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Figure 1: Optical spin orientation induced by transitions across the band gap in a Ga-As semi- 
conductor. 

transition probability of (I) is three times larger than the transition probability of (II). Hence, 
in the undistorted cubic Ga-As, the degree of spin polarization becomes, 



P = - — ^ = 0.5 = 50 % 
3 + 1 



(1) 



The degree of spin polarization can be made 100% with a slight modification in the Ga-As 
material. The spherical symmetry of the crystal field can be destroyed by distorting the cubic 
crystalline symmetry in the electric field of Ga-As through epitaxial growth on a lattice mis- 
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matched substrate, and accordingly, the energy levels get splitted with mj = ±3/2 having a 
higher energy ~ 0.1 eV than nij = ±1/2 (see Fig. [1]). By tuning the photon energy to the 
edge of the band gap, fully polarized spin up or spin down electrons can be produced in the 
conduction band by simply switching between LCP and RCP lights. 

Here it is important to note that, the surface of the Ga-As material can be made with negative 
electron affinity. It manifests that the electrons at the bottom of the conduction band can escape 
into vacuum without any additional supply of energy. 



3 Formal description of spin polarized electrons 
3.1 Quantum behavior of a single spin 

According to the fundamental postulate of quantum mechanics, a wave function ^(r, asso- 
ciated with a microscopic particle having linear momentum p and energy E can be expressed 
as, 

V'(f,t) = A e^^^-^^-"^*) (2) 

where A being a constant. 

uj and k satisfy the following relations, 

nLO = E (3) 

and, 

nk=p (4) 

The phase velocity of the wave packet is u/k. 

Since we are dealing with non-relativistic electrons, the space and spin parts of the wave 
function can be written separately. Thus we can write, 

^l;{f,t,s) = xis)Ae'('^-'^~'^'^ (5) 

Due to the presence of a magnetic field or an effective magnetic field, if there exists an energy 
difference between the two spin states of an electron i.e., AE {= E^ — E^) ^ 0, then the two 
spin states will travel with distinct phase velocities. 

In the forthcoming two subsections (3.2 and 3.3), we discuss the behavior of a single spin 
polarized electron and an ensemble of spin polarized electrons, the so-called a spin polarized 
electron beam. A spin polarization vector P is introduced for the description of a spin polarized 
beam, where its component for the quantization direction ^ is determined by the number of 
up (t) and down (|) spins along 5^. For the complete specification of a spin polarized electron 
beam, the polarization vector P as well as the number of electrons has to be specified. 

To be precise, the main goal of studying a spin polarized electron beam is to answer a funda- 
mental question, which is, "If a beam of n spin polarized electrons described by a polarization 



6 



vector P, traverses through a sample of magnetization M which is oriented at an arbitrary angle 
with respect to P, what will be the magnitude and direction of P, when the electron beam will 
emerge from the sample ?" 



3.2 Single electron polarization in Pauli spin formulation 

Let us start with the basis vectors in cj^ diagonal representation, viz, 

= I T> = J ^ (6) 

Xdown ^ I ^ 1 ^ 

For an electron with spin in arbitrary direction, the spin part of the wave function can be written 
as, 

-(:) 

^ ^ " ( ) ^ ( 1 ) 

where, ui and U2 being arbitrary complex numbers. 
By definition, 

X* = {ul u*^) (10) 
Using the normalization condition x*X = 1) we get, 

U^Ul + U2U2 = 1 (11) 

This is just the one constraint, there are still three free choices of parameters for the determi- 
nation of ui and U2, since these are the complex numbers. 

Expectation values of spin components along x, y and z axes are given by {ax), {(Jy) and (cxz), 
respectively, where the components are. 

These are the so-called Pauli spin matrices. Now considering the general spin wave function Xi 
the expectation values of the Pauli matrices are as follow: 

(o-x) = ixWxlx) 

= ulu2 + U2U1 (13) 
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{(yy) = {X\(yy\x) 

= i[u\ui — U[U2\ 



(14) 



= u\ui — U2U2 

The polarization vector P of a single spin is defined as, 



P 



(15) 
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(16) 



Therefore, using the Eqs. (jl3p . (|14p and (jlSp . the polarization vector P can be written as. 



P 



( pA 


( 


Py -= 




\ P^ J 





( [u{U2 + U*2Ui) \ 
i{u2Ui — ulU2) 



(17) 



Thus, depending on the choices of Xj the polarization P will have some specific forms. For 
illustrative purposes here we set few examples. 

Case 1: The spin is directed along +X axis. 
For this particular choice, we can write. 



1 / 1 



X 



(18) 



V2\l J 

Using this expression of X; the expectation values of the Pauli spin matrices are, 

{a,) = l, {ay) = 0, {a,)=0 (19) 
So the polarization vector P becomes, 

( P.\ (l\ 

(20) 
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Case 2: The spin is directed along —X axis. 

For this specific case the wave function can be chosen as, 



X 



72 I -1 



(21) 



With this form of X) we get, 

{a^) = -1, {ay) = 0, {a,) = 
Therefore, the polarization vector P becomes. 



(22) 









p = 


Py - 









V ) 



(23) 



In a similar way, we also get the polarization vector (P) when the spin is aligned in other 
directions. For representative examples, the results are summarized in Table [H 



Table 1: Spin wave function and polarization vector for spin alignments along X, Y and Z axes. 



Spin alignment 


Spin wave function 


Polarization vector 


Along + X axis 
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' -1 ^ 
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Along + Y axis 
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(;) 
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Along — Y axis 
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' ^ 
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Along + Z axis 


( 
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Along — Z axis 
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' ^ 


) 



Determination of the complex coefficients ui and U2 



9 



The complex coefficients ui and U2 can be derived, in general, if the polarization vector P is 
given. The vector P can have any arbitrary direction in space. The factors and |u2p 

correspond to the probabilities of getting the spin along +Z and —Z axes, respectively, and 
accordingly, they {ui and U2) satisfy the relation + = 1. 



Let us start with the general spin wave function. 



X 



Ui 
U2 



COS a e 
sin a 



i/3 



(24) 



Using this spin wave function, the expectation values of the Pauli spin matrices can be calculated 
those are in the form: 



COS a e *^ sin a 
sin 2a cos /? 



1 

1 



cos a e 
sin a 



if3 



(25) 



COS a e *^ sin a 
— sin 2a sin /? 



-i 

1 



cos a e 
sin a 



(26) 



ixWzlx) 

cose '=~*^ 
COS 2a 



1 

-1 



cos a e 
sino 



if3 



(27) 



Therefore, the polarization vector P is. 



P 



From the above matrix equality we get, 





p. 


\ 




^ sin 2a cos 2/3 \ 




Py 






— sin 2a sin /3 


V 


Pz 


) 




^ cos 2a J 



(28) 



sin 2a cos f3 = Px 
- sin 2a sin (3 = Py 
cos 2a = Pz 



(29) 
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The above three expressions in Eq. ([2^]) can be solved to get the values, 



sma 



cos a 



COS /3 



P. 



I -P. 



(30) 



With these expressions, the components ui and U2 of the spin wave function can be written in 
terms of the polarization vector components as, 



Ul 



I- P. 



U2 



(31) 



Therefore, the probabilities of getting the up (t) and down (|) spins along +Z and —Z axes 



Spin Up Probability 




Figure 2: Coherent decomposition of a single spin along the quantization axis Z. 



become, 



U2 



l + Pz 



I -P. 



(32) 



which obey the relation + |^^2p = 1- The result is schematically described in Fig. [2l 
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3.3 Spin Polarized Electron Beam 



In actual practice we generally deal with an ensemble of spin polarized electrons, instead of a 

single polarized electron, and accordingly, the formalism described above needs to be extended 

for the description of an ensemble. If all the individual electrons are aligned to a particular 

direction, the beam is said to be completely polarized, while for a partially polarized beam the 

individual electrons are aligned in different directions. So in short we can write, 

|P| = 1: pure spin state. 

< |P| < 1: mixed spin state. 

For the complete description of a spin polarized electron beam, both the number of electrons 
and the spin polarization vector (P) have to be specified. Here a basic question that arises is, 
how one can determine P for an ensemble of spin polarized electrons ? 

Now each individual electron i (say), having its spin orientation along ^ axis can be completely 

described by a linear combination of two components n| ^ ^ ^ and ^ ^ along ^ and — ^ 

axes, respectively. Thus, for the entire electron beam, the probabilities of getting the up and 
down spin electrons along ^ are, 

nt = ^|ni|2, ni = ^|4|2 (33) 

i i 

The degree of spin polarization P^ of a spin polarized electron beam relative to a quantization 
axis ^ is defined as, 

P, = ^^1^ (34) 

where, | means along +^ axis and | means along — ^ axis. In terms of intensity of magnetization 
(I) the spin polarization can also be defined which is in the form, 

rmaj jmin 

P = —-. (35) 

The +ve value of |P| manifests that the polarization is along majority spin direction, while, the 
— ve value of \P\ means the polarization is along minority spin direction. 



3.4 An alternative formulation for the description of spin polarization 

The polarization vector P can also be written in the density matrix formulation as. 





( p.\ 


( Tr{a.,p) \ 


p = 


Py - 


= Tr{ayp) 
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where, p is the polarization density matrix which can be expressed in the form, 

.■if ^ + '"'-'''A (37) 



3.5 Spin analyzer 

Spin analyzer is a particular device that is used for the measurement of spin polarization of an 
electron beam. It consists of a spin filter and an electron detector. Only the intensities of 'spin 
up' and 'spin down' electrons relative to the quantization axis can be measured with the aid 
of a spin analyzer. It does not measure their phases. Generally three parameters are used to 
measure the properties of a spin analyzer. They are as follow. 

(i) A direction of maximum transmission, denoted by the unit vector e in space. 

(ii) A transmission factor T^"^ for electrons that are fully polarized along e. 

(iii) A transmission factor T^**^ for electrons that are fully polarized along — e. 

Now to reveal the polarization of incident electron beam with an analyzer, the response of the 
analyzer has to be known. The response is completely described by two factors, Tr and ATr, 
where, 

rp ^ frpmax I rpmin\ 

atr = ^(T^'^'^ -r^^") (38) 

Choosing the quantization direction as the maximum transmission direction, spin analyzer can 
be completely described by the transmission filter matrix as, 

1+ AT„ 



F = Tn\^^^'^- 

( Tr + \ATr 

= 




(39) 



If Tr and ATr are known parameters, the polarization of the incident beam along ^ axis can 
be measured by two successive experiments. Now the incident beam and the transmitted beam 
are related by the expression, 

= Tr{poF) (40) 

-'0 



where. 



Pr,^ + iP^o 1 - PgO 



po = ^( ) (41) 
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Therefore, we can write, 



Tr + -P^^ATjf^ 
2 



/ 

j-Tn 



(42) 



3.6 Spin Filter 



An electron beam becomes spin polarized after passing through a spin filter. The transmitted 

beam becomes polarized even if the incident beam is unpolarized. 

Let the incident and transmitted beams are characterized by, {Iq, pq, P^) and {I,p,P). 

Let us consider an incident unpolarized electron beam with = 0. So, P^ = P^ = P^ = 0. So 

the polarization density matrix becomes. 




(43) 



So, we can write. 



To 



Tr{poF) = -Tr{F) = Tr 



(44) 



Now, the polarization density matrix for the transmitted beam is defined by the relation, 



P=jPoF= — 



So, 



1 + ^ 



Comparing Eqs.[lT]and[36]we can write, 



= = and P^ 
So the polarization matrix becomes. 



•2Tr 



(45) 



(46) 







( ° 


\ 


p = 













V 2Tfl 





(47) 



Therefore, the incident unpolarized beam has got some polarization along the quantization 
direction after passing through the spin filter. Tr and ATr can be calculated by determining 
Tmax and Tmin where they can be evaluated through experiment, by the same way that is used 
to calibrate 'Mott polarization detector'. 



14 



3.7 Fundamental interactions of a spin polarized electron beam with matter 



Experimental observation suggests that when a spin polarized electron beam traverses through 
a nonmagnetic noble metal, like Cu, the polarization of the beam becomes reduced, while the 
spin polarization gets increased when the spin polarized electron beam traverses through a 
ferromagnetic metal. It is generally suggested that these changes are due to 'spin flips', i.e., 
transition from up spin state to down spin state and vice versa. In actual practice, the change in 
spin polarization is not typically due to spin flips, but there are some other processes which are 
responsible for this mechanism. For illustrative purposes here we mention some of them as follow. 

• Spin precession: It is associated with a phase shift between the up and down spin 
wave functions which provides a precession of P about the quantization axis, keeping the 
magnitude (|-P|) as invariant. 

• Spin exchange scattering: This is related to the transition from one spin channel to the other 
spin channel. The polarization vector P rotates relative to the quantization axis ^, with a 
change in the component P^. 

• Spin dephasing: It is related to the out of phase precession of the individual spins about 
the quantization axis, and accordingly, the magnitude of the component P± (perpendicular to 
the quantization direction) gets decreased. In non-magnetic materials, this effect may lead to 
complete spin randomization. 

• Spin selective absorption: This particular phenomenon originates from the preferential 
scattering of spin in one spin channel. It causes a rotation of P relative to the quantization 
axis resulting a change in magnitude, i.e., |P| 7^ l-P*^!, where jP*^! is the initial amplitude of the 
polarization vector. 

3.8 Two Current Model 

The total number of electrons in each atom those are responsible for electrical conduction, 
n{EF), is directly proportional to the density of states at Ep. This proposes that Cu to be 
a bad conductor, while, Fe, Co and Ni having a large value of n{Ep) to be good conductors. 
But in actual practice, Cu has the highest conductivity among all of these transition metals. 
Naturally, a question arises why there is such a discrepancy ? 

This discrepancy was removed by Neville Mott. He guessed that 's' electrons carry electric 
current, and, due to scattering of electrons from 's' band to 'd' band, the electrical resistance 
appears. Therefore, stronger scattering takes place when more 'd' states are available. Since 
this scattering is 'spin selective', only the electrons having the same spin can jump from 's' band 
to 'd' band. Spin mixing between the two channels are forbidden. 

In the 'two current model', the high resistivity of the transition metals with a partially filled 
'd' shell is thus produced by scattering of 's' electrons to the empty 'd' states. The modern band 
theory supports the strong hybridization of 's' and 'd' electrons, and accordingly, the excitations 
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from 's' to 'd' states arise naturally. 

In the 'two current model', the two spin states carry the electric current separately and 
independently of each other. This is due to the spin selective scattering process of 's' electrons 
to the vacant 'd' band, which is the key aspect of this model. The basic feature of the 'two 
current model' is based on 'Fermi's Golden rule' which states that the scattering probability of 
the conduction electrons, that leads to an increase in resistivity, is proportional to the density 
of the vacant 'd' states available above the Fermi energy. 

Since the conduction in two different spin channels occurs independently, the net conductivity 
becomes, 

a = a^+ai (48) 



where, cj-f 



and a 



""-^^ '^^ . This 'two current model' is valid in diffusive as well as in 



Majority 
Spin 




d band 



s band 



Spin Excitation in the "two current model" 
Figure 3: Spin excitations in the two current model. 

ballistic regime. 

So, in a nutshell what we can say about 'two current model' is: 



The two current model predicts that, electric current flows in two separate spin chan- 
nels independently, and the total conductivity is obtained by adding the individual 
contributions of the two separate spin channels. The resistivity gets increased when 
the electrons in one of the spin channels undergo spin conserving scattering from 's' 
band to the empty states of 'd' band. 
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The 'two current model' is illustrated schematically in Fig. O 



4 Formal theory of spin transport 

4.1 Spin accumulation voltage 

Let us consider a steady current is flowing along the direction x through the interface between 
a ferromagnetic metal (FM) and a non-magnetic metal (NM). The quantitative description of 
this current can be formulated in terms of 'spin dependent chemical potential'. 



The 'spin dependent chemical potential' of majority electrons of the system is specified by 
the energy change which is associated with the addition of one majority electron to the 
system. 



If no current flows through the interface, Fermi levels of both the metals (FM and NM) will be 
rearranged to the identical value /io, where, 

MO = (49) 

According to Ohm's law, when a current flows from FM to NM along x direction, we can write, 

= --3 (50) 
ox a 

where, e is the electronic charge and a is the total electrical conductivity due to up and down 
spin electrons, i.e., a = + a^. Spin flips, required to make transitions from one conductivity 
channel to the other, can be neglected at low temperatures. But, this is a good approximation 
even at room temperature, since the spin flip length is larger than the electronic mean free path 
at that temperature. 

Now interesting spin dependent phenomena must occur at the interface between a ferromag- 
netic metal, characterized by a preferential conduction in one spin channel, and a non-magnetic 
metal having equal conduction for both the spin channels. From the charge neutrality condition 
at the interface, we get, j = constant. Therefore, in fi-x curve, /ig becomes a straight line with 
a change in slope at the interface, as cjat ^ ap, where ajy and ap correspond to the conduc- 
tivities of the non-magnetic and ferromagnetic metals, respectively. But the description is not 
complete because of the fact that, in the ferromagnetic metal unequal number of up and down 
spin electrons provides different electrical conductivities associated with each spin channel. 
Let us define, 

T 

CTp = (1 - ap)ap (51) 



17 



where, ap is a parameter lying between to 1. Now, due to the asymmetry of conduction in 
the two spin channels, ap 0.5. In a bulk ferromagnet, minority spin channel has the lower 
conductivity, and therefore, ap > 0.5. 

Asymmetry in conduction through the two spin channels results asymmetric current. Now we 
define, 



JF 



JF-" 



T 



Pfj 

(1 - PF)j 



(52) 



and. 



JN^ = (1 - Mj 



(53) 



where, Pf and are the two parameters. The spin polarization P of the current is defined as, 

P 



(54) 



where, a corresponds to F{N) when the spin polarization is defined in case of FM(NM). When 
the spin polarized current enters from FM {I3p / 0.5) to NM {(3j\f = 0.5), / 0.5 and P / at 
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Figure 4: Variation of spin dependent chemical potentials across the interface, when a steady 
current is flowing from a ferromagnetic metal to a non-magnetic metal. 



the interface. This is due to the fact that, spin flips are rare, so it will take time to equilibrate to 
the value Pn = 0.5 away from the interface. This introduces a voltage drop at the interface, the 
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so-called 'spin accumulation voltage' (Vas) as shown in Figs. [Hand O The chemical potentials 
due to up (I) and down (|) spins, fi^ and fi^, must be continuous across the interface. For 
a non-magnetic metal since the two spin channels conduct equally, /Uq will be the average of 
fi^ and fi^. On the other hand, it is a weighted average for a ferromagnetic metal. Thus, for 
non- magnetic metals (NM), 

(55) 



1 



and for ferromagnetic metals (FM), 



fJ'OF = CtpfJ'^ + (1 - OiF)f^^ 



(56) 



When current flows from FM to NM, V^s can be written as. 
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Figure 5: Variation of spin dependent chemical potentials across the interface, when a steady 
current is flowing from a non-magnetic metal to a ferromagnetic metal. 



cVas = ^J'OF{x = 0) - hon{x = 0) 

1" 



(57) 



In this particular case, ap > 0.5. Thus, Vas / and /i^(0) > At^(O). For the reverse situation, 
i.e., when the current flows from NM to FM, we can write. 



eVAs = fJ'ONix = 0) - fiorix = 0) 



(58) 



Hence, in this situation, fj.^{0) > /x^(0). 
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To summarize, the spin accumulation voltage (Vas) can be defined as the jump in spin 
averaged chemical potentials, developed at the FM/NM interface, because of the asymmetric 
conduction of the two spin channels through the interface. 



4.2 Diffusion equation for spins 

A spin voltage is defined as the difference between spin up and down chemical potentials and it 
reflects the difference in the number of electrons in the two current carrying spin channels. The 
diffusion equation for the up (j) and down (J,) spins can be written in terms of their chemical 
potentials as, 



^ , " ^ = (59) 



Solving Eq. [59] we get. 



T i 



[/iT(0) - /i^(0)]e-t (60) 



where, A is given by, A = \/DTse- The decrease of spin voltage (//^ — iJ,^) with distance from the 
interface is determined by the characteristic equilibrium distance A, the so-called 'spin diffusion 
length'. Thus the spin dependent chemical potentials (/i^ and fi^) equilibrate away from the 
interface (x = 0) over a characteristic length (A) and a characteristic time (r^e). In order to 
observe the spin dependent transport through a conductor, the sample size should be less than 
or comparable to the 'spin diffusion length'. 



5 Spin transport at atomic scale 
5.1 Introduction 

Recent progress in nanoscience and nanotechnology has enabled us to develop various magnetic 
materials and magneto devices which reveal several fascinating phenomena. This new era in the 
study of magnetism has been initiated in 1988 after the discovery of Giant Magneto Resistance 
(GMR) effect [2] observed in magnetic multilayers formed by alternating magnetic and non- 
magnetic materials. 

In the absence of an external magnetic field, the exchange coupling between adjacent mag- 
netic layers through the non-magnetic one aligns the magnetization vector anti-parallel to each 
other. Then, when a magnetic field, strong enough to overcome the anti-ferromagnetic coupling 
is applied, all the magnetization vectors align along the field direction. This new parallel configu- 
ration exhibits an electrical resistance which is much smaller than that of the anti-ferromagnetic 
configuration. This drastic change in resistance is known as the GMR effect which is illustrated 
in Fig. [6l Mathematical definition of GMR is given by. 
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In normal GMR effect, Rap > Rp, and R is unbounded. Another definition used for GMR is, 



R' 



Rap - Rp 
Rap 



(62) 



In this particular case, < i?' < 1. 

Inverse GMR effect has also been found in some materials for which Rp > Rap- The study of 
GMR effect is of fundamental importance as it established the fact that the spin of an electron 
can also take an important role in transport phenomena. Today, the study of GMR effect has 
drawn much attention from academic circuit to commercial levels [3] as GMR based magnetic 
data storage devices are used in every computers. Not only that, spintronics in low dimensional 
systems provides some real advantages over bulk metals and semiconductors. In molecular 
systems, the conventional mechanisms for spin decoherence (spin orbit coupling, scattering from 
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Figure 6: Schematic representation of GMR sensor. 

paramagnetic impurities, etc.) get reduced. Hence, in nanoscale systems, we can expect the 
spin coherence time to be several orders of magnitude larger than in bulk systems. Therefore, 
the study of spin dependent transport and spin dynamics is of great importance to understand 
and to develop the field, spintronics. 

Spin dependent transport through mesoscopic systems has attracted rich attention in the 
development of spintronics [H O El [7] . Several experimental efforts have been made to study 
spin transport in quantum dots (QD) [3 [9l [lOl [11] and molecular systems pT j [T2 | [T3 l [T^. 

In these systems, electrical resistance depends on the spin state of electrons passing through 
the device, and it can be controlled by an applied magnetic field. This is due to the imbalance 
in transmission probability for spin up and down electrons through the device |15l [IB] . 

The spin dependent transport through nanostructures such as a QD system, can be inves- 
tigated by coupling it to magnetic or nonmagnetic leads and passing a current through the 
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system [TH [T71 [TSl [THl [2D]. Several theoretical as well as experimental studies have been made 
on the spin transport through quantum dot devices |21|, [22} [23] . 

Aim of the model that has been analyzed below is to study analytically the coherent spin 
dependent transport of a magnetic QD system attached to nonmagnetic leads, with the magnetic 
quantum dot (MQD) device having an arbitrary magnetic configuration. This analysis has been 
done using Transfer matrix method in the framework of tight-binding Hamiltonian |24t \2U[ 
[271 [28]. 



5.2 Description of the model and formalism 

The model chosen is a heterostructure formed by a sequence of magnetic and nonmagnetic sites 
connected to nonmagnetic leads. The spin dependent electronic transport through this structure 
has been studied including the spin flip scattering effect. The heterostructure is basically a 
repetition of a chosen unit cell formed by magnetic and nonmagnetic atoms. The geometry 
of the NM/MQD/NM structure is shown schematically in Figs. [7| and [51 The direction of 
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Figure 7: A magnetic quantum dot device (magnetic spacer) is attached to source and drain. 
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Figure 8: Schematic view of the model. A magnetic conductor is connected to nonmagnetic 
source and drain. 

magnetization on each magnetic site in the MQD device is chosen to be arbitrary, specified in 
each site n, which in spherical polar coordinate system is defined by two angles On and (pn- 
Here On is the angle between magnetization direction and the Z axis, (pn is the azimuthal angle 
of magnetization measured from X axis at site n as shown in Fig. [9l The whole system can be 
described by a Hamiltonian of the general form, 

H = Hl + Hld + Hd + Hdr + Hr (63) 
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The spin polarized electrons for a A^-site quantum dot can be described within the effective 
one-electron approximation in tight-binding framework and the nearest-neighbor approximation 
as, 

Hd = ^cl (^en - hn.a) C^ + Y^ (cl,tn,n+lCn+l + c|^+itn,n+lCn) (64) 



where, d = ( 4^ 



, / tn,n+l \ 

tn,n+l ~ n 

V tn,n+l ) 

The term e„ of Eq. [M] describes the on-site energy. The term hn-S describes the interaction 
of electrons with the magnetic atoms in MQD device. This term allows the spin flips at the 
magnetic sites. /i„ is the amplitude of the spin flip parameter at site n. a is the Pauli spin 



Z A 




operator having components {ax,cry,crz)- Spin flip scattering is dependent on the magnetic 
moment orientation of the atoms in MQD device with respect to the Z axis. The term tn,n+i 
is the nearest-neighbor hopping integral, c^na{cna) is the creation(annihilation) operator for an 
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electron on site n with spin a. The Hamiltonian of the left(right) lead is defined as, 



Hl{R) = J2 ^^iCi + J2 (citL(R)Ci+l + c|'_^itL(R)Ci) (65) 



where tL{R) refers to the hopping integral between the sites of the left (right) lead. Finally, the 
Hamiltonian that corresponds to the coupling of the MQD device to the leads can be expressed 
in the form, 

HlD{RD) = (co(N_,_i)tLD(RD)Ci(N) + c|(-]v^)tL£,(j^D)Co(N+l)) (66) 

Throughout this study, it has been assumed that the nonmagnetic leads are ideal, i.e., their 
resistances have been neglected. The main contribution to the resistance in this ballistic device 
has come from contact resistance and spin scattering in MQD. 

Using this model the spin dependent transmission coefficients are calculated to investigate the 
transport properties, in transfer matrix formalism. In this calculation the Schrodinger equation 
is written in terms of ipn,^ (localized Wannier basis), which express the amplitude of the wave 
function at site n, with energy E and spin a. 

5.2.1 Transfer Matrix method 

In order to calculate the spin dependent transmission probabilities, the Schrodinger equation 
for the electron wave function must be solved in the MQD device. So the starting point is the 
time independent Schrodinger equation in the MQD device, which can be written as, 

H\ct>) = E\<^) (67) 

where, 

!</') = E[^*Tl^T) + ^alU)] (68) 

i 

Here, \(f)) is written as a linear combination of spin up and spin down Wannier states. Now, 

h.a = hxCx + hyay + /i^cr^ 

/ h cos On h sin One""^^" \ 
\ h sin Onc'^'^" -h cos On J 

So the Hamiltonian for the whole MQD device can be written as, 
^ _ / t A (^n]-hcosOn -/i siu 6'„e~*'^" \ / CnT \ I 

t ,w * ^ \ ( ^^"+1'^ V r.t \( ^ ^ \ ( C".T 



^ V t J \ Cn+l,l J ^ ^ \ t J \ Cn,l 



E 

n 

Hi + H2 (70) 
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where, 



Hi = X! I" t) (cnT - /icos6'„) (n t I t)/isin6'„e '"^"(n i 

n n 

DhsinOne'^'^in] | + ^ |n j) (e„| + /icos6'„) (n j | 

n n 

H2 = ^|?iT)t(?i + i,T | + + 1 + 

n n 



(71) 



(72) 



Operating H on \(j)) we get the following two equations relating the Wannier amplitudes on site 
n of the MQD device with the neighboring nil sites, 



{E - e.„ + hn cos 0n) tpn^ + K siu One "^"tpnl = tlpn+1,1 + tlpn-1,1 

hn sin 6'„e*''''" VnT + (E - e„ - /i„ cos 9n) ipni = ^V'n+ia + ti'n-1,1 



(73) 
(74) 



Transfer matrix for the nth site relates the wave amplitudes of nth site with that of (n — l)th 
and (n + l)th sites. So we can write the transfer matrix equation for the nth site as. 



So the total transfer matrix for the whole MQD device can be written as, P = Yl Pi 

l=N 

where, 
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(75) 



Pn 



t t 

t t 



-1 
-1 




(76) 



1 
\ 1 y 

For the whole lead-MQD-lead system, the transfer matrix (T) that relates the wave amplitude 
from the left lead ((-l)th and 0th sites) to the right lead {{N + l)th and {N + 2)th sites) is 
given by. 



'4'N+2l 

V V'iv+u J 



V V'-u J 



(77) 



where. 



T = Mr.P.Ml 



(78) 
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Here, Mji = the transfer matrix for right boundary, i.e., (N + l)th site. It relates the wave 
ampUtudes from A^th site (in the MQD device) to {N + l)th site (in the right lead). Similarly, 
Mx, = the transfer matrix for left boundary, i.e., 0th site. It relates the wave amplitudes from 
1st site (in the MQD device) to (— l)th site (in the left lead). 



5.2.2 To calculate Ml and Mr 

Let us calculate Ml first. For the left (or right) lead, h.a = 0. So the Hamiltonian for the lead 
(left or right) can be simplified as, -ffiead = + where, 

H, = J2i\n])er,^{n]\ + \nl)enia\) (79) 

n 

and 

H2 = Y. (1^ T)i(n + 1, T I + |n l)t{n + 1, i | + |n + 1, ])t{n j | + |n + 1, l)t{n [ \) (80) 

n 

We set en = eo for all n in the leads. Now, from the time independent Schrodinger equation, 
operating H o'n.\(f)) we get two equations for 0th site, relating the wave amplitudes of 1st and 
— 1th sites, as given below. 

[E - eo) -00? = tLoi^i^ + tLip-i] 

(E-eo)0Oi = tLDi^n + tL^P-U (81) 
Now in the lead, according to the tight-binding model, 

iPn = Ae'^""" (82) 

So ip-i can be written in terms of V'o as, 

V'-i = V'oe-*^^ (83) 

where, Pl = ka and 

E = €o + 2tL(R)Cosf3L(R) (84) 
From Eqs. [81] and [82} ipi can be expressed in terms of ipQ as, 

V'l = (^e'^') V'o (85) 
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Hence, the transfer matrix matrix equation for the 0th site is of the form, 

^ V'lT \ ( l^^e^''^ W V^oT \ 

i^n ^ J^e^^^ V^-oi 

V^OT e*^i -f/'-n 

\M J \ e'f^^ J \ Vj_u J 

Similarly, from the time independent Schrodinger equation, operating H on \(p) we get two 
equations for (N + l)th site, relating the wave amplitudes of A^th and (N + 2)th sites as given 
below. 



{E - eo)V'Af+iT 
{E - eo)V'Ar+U 



(87) 



Exactly as before the transfer matrix matrix equation for the {N + l)th site can be written as. 
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(89) 



(90) 



5.2.3 To calculate the transmission probability of up and down spin electrons 

Thus the transfer matrix for the whole system (lead-MQD-lead) can be written as, 

T = matrix for (N+l)th site. product matrix for N-th site to 1st site, matrix for 0-th site 
= Mr.P.Ml (91) 
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In order to calculate the transmission coefficient for the incident electrons with up or down spin 
using the transfer matrix method, the wave function amplitudes (Wannier amplitudes) have to 
be specified on proper atomic sites. 



Case 1: Up spin incidence in the left lead. 



The eigenvalue equation involving the transfer matrix relating the Wannier amplitudes from 

sites {N + 2) and {N + 1) to sites and -1 is given by Eq. [771 

Let, 

p^^ = Reflection amplitude for up spin (|) reflected as up spin (|). 
p^^ = Reflection amplitude for up spin (t) reflected as down spin (|). 



-Ti 



Transmission amplitude for up spin (t) transmitted as up spin (|). 
Transmission amplitude for up spin (t) transmitted as down spin (|). 



Now, for the left lead, the wave amplitudes at the sites and -1 can be written as, 



■001 
0oi 



(92) 



Similarly, for the right lead, the wave amplitudes at the sites {N + 1) and {N + 2) become, 



'>PN+21 



(93) 



Therefore, the transfer matrix eigenvalue equation can be rewritten in terms of these wave 
function amplitudes as. 



f ^TTe«(^+2)/3fl \ 

-npi{N+l)f3R 



\ 



( 



1 + pTT 
pU 



\ 



\ 



(94) 



Solving the above equation we can get the values of r^^ and t^^ . The transmission coefficients 
T^^ and T^^ are defined as the ratio of the transmitted flux to the incident flux as, 



-TT 



tR sin/?/; 
Il sin/3i 



-TT|2 
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tL sin ' 

Therefore, the total transmission probabihty for spin up is, 



(95) 



(96) 



Case 2: Down spin incidence in the left lead. 



Similarly, for a down spin incidence in the left lead, the amplitudes at sites —1, 0, (A^ + 1) and 
(A^ + 2) are given by. 



1 + 



(97) 



and 



V'Af+2T 
V'Af+lT 



^ll^i{N+l)betan 



(98) 



So, as before, the transfer matrix equation for down spin incidence can be rewritten in terms of 
these wave amplitudes as. 



( 



V 



1 + 



\ 



j^pllf^ifiL J 



(99) 



This equation can be solved to compute t^^ and r-*^!^. 

The transmission coefficients T^^ and T^'^ are defined by the ratio of the transmitted flux to 
the incident flux as, 



Hi 



*Rsin^|^^|2 

tL sin (3l 

tR sinPn 2 



(100) 
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Therefore, the total transmission probabihty for spin down is, 

n=ru+nT (101) 

At low temperatures, the transport of such mesoscopic systems is ballistic and coherent. In 
this regime, transmitted spin dependent current flowing through the interacting region with 
arbitrary magnetic configurations, symmetrically attached to the two NM leads is given by the 
relation [Ml ED], 

+00 

laa' (^) = ^ / ifs - Id) T^a' (E) dE (102) 

—00 

where, fs{D) = / ~ /^5(D)) gives the Fermi distribution function with the electrochemical 
potential ^s(D) = Ep ± eV/2. Thus, /o-o-' is calculated by determining Tg-o-' from Eqs. [95] and 

[Ml 

5.3 Numerical results and discussion 

In order to investigate the spin transport properties through the MQD device in NM/MQD/NM 
heterostructures, we perform some numerical calculations. To do this we assume that the MQD 
device consists of several unit cells (see Figs. [TOl [T3l [TBI and [T9|) that are made up of some 
magnetic (A) and nonmagnetic (B) atoms having different magnetic moment orientation. By 
applying an external magnetic field, the magnetization direction of the NM/MQD/NM system 
can be switched. 

To illustrate the results, let us first mention the values of the different parameters those are 
considered throughout our numerical calculations. For simplicity, all the on-site energies both in 
the two leads and in MQD device are taken to be zero. The nearest-neighbor hopping integral t 
in the MQD device is set to 1, and, in the two side attached leads the hopping strength tL(_R) is 
also fixed at 1. The Fermi energy Ep is fixed at 0. The spin fiip parameter is chosen as h = 0.5. 
To narrate the coupling effect [311 [32l [331 EH , throughout the study we focus our results for the 
two limiting cases depending on the strength of the coupling of the MQD device to the source 
and drain. Case A: The weak-coupling limit. It is described by the condition t/^(_r) « t. For 
this regime we choose tl = = 0.2. Case B: The strong-coupling limit. This is specified 
by the condition tj;^(/?) ~ t. In this particular regime, we set the values of the parameters as 
tl = Tfj = 0.8. For the sake of simplicity we choose the unit c = e = h = 1. 

Case 1: 

Let us start with a periodic MQD system in the parallel magnetic configuration. The MQD de- 
vice is made up of several unit cells, where each unit cell contains four magnetic (A) atoms, each 
having parallel magnetic moment orientation. The schematic view of the unit cell configuration 
is given in Fig. [TOl Here we study the variation of the conductance 5 as a function of the injecting 
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Figure 10: Unit cell configuration 1 made up of four magnetic atoms A, each having parallel 
magnetic moment. 



electron energy E, both for the weak- and strong-coupling limits, respectively. The results are 
focused for the two separate cases, respectively. One is the transmission of up spin (|) as up 
spin (I) and the other is the transmission of up spin (|) as down spin (|), i.e., the case of spin 
flip. The nature of the conductance spectra for a typical MQD device with total number of site 
= 64 in the limit of weak-coupling is shown in Fig. Ill[ Figure [TTTa) predicts the conductance 
due to the transmission of up spin as up spin. On the other hand, the conductance involving 
only the spin flip effect, i.e., transmission of up spin as down spin, is given in Fig. [TlTb). From 
Fig. IllTa) it is observed that, for some particular energies the conductance exhibits fine resonant 
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Figure 11: g-E spectra for a MQD device with the unit cell configuration 1, in the weak-coupling 
limit, (a) g^i and (b) g-^^. 

peaks, while the conductance drops to zero for other energies. At the resonant energies where 
the conductance approaches the value 2, the transmission probability T goes to unity, since the 
relation g = 2T follows from the Landauer conductance formula, g = 2{e^ /h)T, with e = h = 1. 
All these resonant peaks are associated with the energy eigenvalues of the MQD device, and 
accordingly, we can say that the conductance spectrum manifests itself the electronic structure 
of the MQD device. Now, Fig. [TlTb) shows zero transmission probability for transmission of up 
spin as down spin, for this type of unit cell configuration. The effect of coupling is quite interest- 
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Figure 12: g-E spectra for a MQD device with the unit cell configuration 1, in the strong-coupling 
limit, (a) g^i and (b) 



ing. For our illustrative purposes, in Fig. [12] we plot the g-E spectra for the same MQD device 
{N = 64) in the limit of strong-coupling, where (a) and (b) correspond to the same meaning as 
in Fig. [TTJ Figure [T2ra) predicts that, all these resonant peaks get substantial widths compared 
to the weak-coupling limit. The contribution for the broadening of the resonant peaks appears 
from the broadening of the energy levels of the MQD device in this strong-coupling limit [35] . 
Hence by tuning the coupling strength, we can get the electron transmission across the MQD 
device for the wider range of energies and it provides an important signature in the study of 
current-voltage (I-V) characteristics. Similar to the weak-coupling limit, for this strong-coupling 
case also, there is no spin flipping, and accordingly, we get zero transmission probability for the 
transmission of up spin as down spin. 

Case 2: 

Now we focus on the conductance behavior for a periodic MQD system having anti-parallel 
magnetic moment orientation. The MQD device is made up of several unit cells, where each 
unit cell contains four magnetic (A) atoms, arranged antiferromagnetically. The schematic view 
of the unit cell configuration is illustrated in Fig [T3l For representative examples, in Fig. [TH 



Figure 13: Unit cell configuration 2 made up of four magnetic atoms A, having anti-parallel 
magnetic moment. 
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we plot the conductance-energy spectra for a typical MQD device considering = 64, in the 
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Figure 14: g-E spectra for a MQD device with the unit cell configuration 2, in the weak-coupling 
limit, (a) g^i and (b) g-^^. 
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Figure 15: g-E spectra for a MQD device with the unit cell configuration 2, in the strong-coupling 
limit, (a) g^i and (b) g-^^. 



limit of weak-coupling, where (a) and (b) represent the identical meaning as in Fig. [TTJ Similar 
to the above case, here we also see that for some particular energies the non-zero value of the 
transmission probability of up spin as up spin is achieved. While, the transmission probability 
of up spin as down spin becomes exactly zero for the entire energy range (Fig. [T^ b)). since 
there is no spin flip effect in this cell configuration. On the other hand, with the increase of 
the MQD-lead coupling strength, all the resonant peaks get broadened. The results for the 
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strong-coupling limit are shown in Fig. [TSl where (a) and (b) represent the identical meaning 
as above. Prom both these spectra (Fig. [TUfa) and Fig. [TST a)). we get additional one feature, 
which is the appearance of a global energy gap in the conductance spectrum. The gap appears 
symmetrically about the energy E = with width 2h (since hn = h for all the sites). This reveals 
that the electron conduction takes place beyond the energy h, and no electron conduction takes 
place when the energy E < h. Thus by controlling the strength of /i, we can tune the electron 
conduction. This behavior can be much more clearly explained from our study of current-voltage 
characteristics which we will describe in the forthcoming section. 

Case 3: 

Next, we concentrate our study on the behavior of electron conduction for a periodic MQD 
device which is made up of several unit cells, where each unit cell contains four magnetic (A) 
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Figure 16: Unit cell configuration 3 made up of four magnetic and nonmagnetic atoms A and 
B, A having antiparallel magnetic moment. 
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Figure 17: g-E spectra for a MQD device with the unit cell configuration 3, in the weak-coupling 
limit, (a) g-^-^ and (b) gi^\^. 

and nonmagnetic (B) atoms, arranged alternately. The unit cell configuration is schematically 
given in Fig. [TBI The variation of the conductance g with the energy E for a typical MQD 
device with = 64 in the limit of weak-coupling is shown in Fig. [T71 While, the results for the 
strong-coupling case are plotted in Fig. [THJ For both these two cases, (a) and (b) correspond to 
the identical meaning as described earlier. The properties of the resonant peaks both for the two 
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Figure 18: g-E spectra for a MQD device with the unit cell configuration 3, in the strong-coupling 
limit, (a) g^i and (b) gn- 

coupling cases, i.e., sharp peaks in the weak-coupling limit and broadened peaks in the limit of 
strong-coupling, remain unchanged like the previous unit cell configurations. Also, the behavior 
of getting up spin as down spin through this MQD device remains invariant as before. The 
argument of this vanishing transmission probability is exactly similar as we mentioned earlier. 
For this particular cell configuration we get additional one interesting feature which is related to 
the existence of the global energy gaps at the two distinct regions of the energy E (see Fig. [TTTa) 
and Fig. [TSlf a) ) . The widths of these two gaps are exactly identical to each other, which is given 
by h. The existence of a global energy gap manifests that the system acts as an insulator along 
the perpendicular direction in the long length of MQD device. Here we get two such global 
energy gaps, and accordingly, this particular MQD device can be identified as double gap spin 
valve system, that can be used as a mutual spin switch. 

Case 4: 

Finally, we consider the particular case where each magnetic moment is oriented at an angle 9n 
with respect to the vertical Z direction. For convenience, we consider a single unit cell which 
contains six magnetic moments. The unit cell configuration is schematically presented in Fig.ll9[ 
For illustrative examples, in Fig. [20] we plot the conductance-energy spectra for a MQD device 




Figure 19: Unit cell configuration 4 made up of six magnetic atoms with their magnetic moments 
oriented at an angle 9 with respect to the Z axis. 
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Figure 20: g-E spectra for a MQD device with the unit cell configuration 4, in the strong-coupling 
limit, (a) g^i and (b) gn- 6 is chosen as it/3. 




Figure 21: g^i-O (red curve) and g^i-O (blue curve) for a MQD device with the unit cell config 
uration 4, in the weak-coupling limit. The energy E is fixed at 0.25. 



considering = 64 in the limit of strong-coupling. Figure [2CT a) represents the g-E spectrum 
associated with the transmission of up spin as up spin, while Fig. [20jb) corresponds to the g-E 
curve associated with the spin flip effect, i.e., the transmission of up spin as down spin. For 
this unit cell configuration, we assume that all the magnetic moments are aligned at an angle 
9 = tt/S [On = 9 for all n) with respect to the preferred Z direction. Quite interestingly we see 
that, the MQD device with this particular unit cell configuration exhibits spin flip transmission. 
The reason for non-zero spin flip transmission probability is as follow. The interaction of spin 
of the injected electron with the local magnetic moments associated with each site of the MQD 
device is given by the term /i„.(T in the Hamiltonian (Eg. I64p . Unlike for unit cell configurations 
1, 2 and 3, where the moments are taken to be either parallel or anti-parallel to Z axis, in 
this case all the local moments are aligned at an angle 9 with respect to Z axis, i.e., they have 
non-vanishing x and y components. Therefore, the interaction term contains ax and ay and so 
also (74. and cr„, which are responsible for spin flip. Now, to incorporate the dependence of 9 
on conductance, in Fig. [2T]we present the results for a MQD device considering N = 10, where 
the red curve corresponds to the conductance (gn) due to pure spin transmission and the blue 
curve represents the conductance {g-^i) due to spin flip transmission. 
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5.3.1 Current- Voltage characteristics 

All the basic features of spin dependent transmission become much more clearly visible from 
our study of current-voltage (I-V) characteristics. The current through the MQD device can be 
computed by integrating the transmission function T, as prescribed in Eg. 11021 The transmission 
function varies exactly similar to that of the conductance spectrum, differ only in magnitude by 
the factor 2 since the relation g = 2T holds from the Landauer conductance formula. Here we 
describe the I-V characteristics for all these four different unit cell configurations one by one. 

Case 1: 

First, we describe the results for the unit cell configuration 1. As illustrative examples, in Fig. [22] 
we show the variation of the current (/) as a function of the applied bias voltage (V) for a typical 
MQD device with = 8, where the first and second rows correspond to the results for the weak- 
and strong-coupling limits, respectively. In the limit of weak-coupling, current shows step- like 

(b) -------- 





-4 -2 2 4 
V 

(d) -------- 





-4 -2 2 4 
V 

Figure 22: I-V curves for a MQD device with the unit cell configuration 1. The 1st and 2nd 
rows correspond to the weak- and strong-coupling limits, respectively, (a) I-^i, (b) I-^^, (c) I^i 
and (d) I^^. 

behavior as a function of bias voltage (see Fig. [22lfa)). This is due to the presence of sharp 
resonant peaks in the conductance spectrum (Fig. [TO a) ) . as the current is obtained from the 
integration procedure of the function T. With the increase of the bias voltage, the difference in 
the electrochemical potentials of the leads increases, allowing more energy states accessible to 




-4 -2 2 4 
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the injected electron to get passed through the MQD device from the source to drain. Therefore, 
the current increases providing a jump in the I-V characteristics. On the other hand, the current 
varies almost continuously when the coupling strength is increased (Fig. [2^c)). because in this 
coupling-limit the conductance peaks get broadened which provides larger current amplitude. 
Both in these two coupling limits, the current drops to zero for the entire bias voltage due to 
spin flip transmission. 

Case 2: 

Now we consider the case for the unit cell configuration 2. The I-V spectra for the MQD device 



0.4 

(a) 
- 



-0.4 



6 

(c) 




-6 



-4 -2 



-4 -2 




V 




V 



6 
(b) 




-6 



6 

(d) 




-6 



-4 -2 




V 



-4 -2 




V 



Figure 23: I-V curves for a MQD device with the unit cell configuration 2. The 1st and 2nd 
rows correspond to the weak- and strong-coupling limits, respectively, (a) Iii, (b) /fj, (c) /|| 
and (d) I^. 



{N = 8) are shown in Fig. [23l where the upper and lower rows correspond to the identical 
meaning as in Fig. [22j All the basic features due to the coupling are the same as discussed 
in previous case. But there is one additional important feature due to this type of unit cell 
configuration. Because of the presence of an energy gap in the conductance spectrum, the 
non-zero value of current is obtained after overcoming a threshold bias voltage. 

Case 3: 



Here we study the results for the unit cell configuration 3. The characteristic features of the 
I-V spectra for a MQD device with = 8 is illustrated in Fig. [2^ All the essential properties 
are almost invariant as presented in Fig. [22l In this particular cell configuration, the presence 
of the double energy gap in the conductance spectrum ensures the saturation like behavior in 
the current-voltage curve. 
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Figure 24: I-V curves for a MQD device with the unit cell configuration 3. The 1st and 2nd 
rows correspond to the weak- and strong-coupling limits, respectively, (a) I-^-^, (b) I-^^, (c) /n 
and (d) 



Case 4: 

Finally, we come to the case where the unit cell of the MQD device is of the type configuration 




Figure 25: I-V curves for a MQD device with the unit cell configuration 4. Here we set 6 = tt/3. 
The 1st and 2nd rows correspond to the weak- and strong-coupling limits, respectively, (a) 
(b) ^Ti' (c) -^TT ™d (d) 



4. The results are shown in Fig. [25] for a MQD device considering = 8. The new feature 
associated with the I-V characteristics in this case is the non-zero value of current due to spin 
flip transmission. Another important feature is that the magnitude of current due to spin flip 
transmission is much smaller than that of pure spin transmission. 
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6 Concluding remarks 



In this dissertation, we have addressed some issues of spin polarized transport through matter. 
After introducing the new field-Spintronics, in the first section we have discussed about the 
generation of spin polarized electron source and its detailed mechanism. Then we have intro- 
duced the formal theory needed to study the spin polarized transport. Then we have talked 
about the 'Two Current Model' which is the stepping stone towards the development of modern 
spintronics. This model illustrates the fact that the two spin states of electrons conduct current 
separately and independently to each other, giving rise to different resistivities associated with 
each spin channel. 

Next, we have tried to analyze the spin dependent transport in nanoscale systems. We have 
discussed spin transport in different types of NM/MQD/NM heterostructures, including the spin 
flip effect, within the tight-binding formulation. 

The MQD spacer has been considered to be made up of an array of quantum dots, which 
is basically a repetition of unit cells chosen differently in different cases. For the first three 
cases the unit cells consist of magnetic QDs only, whereas in the last case it was thought to be 
periodically juxtaposed magnetic and nonmagnetic QDs. 

In the theoretical formulation we have obtained the general formulae for calculating the spin 
dependent transmission coefficients, which depends on the arbitrary magnetic configuration of 
the MQD device. 

Depending on the local magnetic moment orientation on each site, in the MQD device we 
obtained gapless, single-gap, double-gap conductance spectra, respectively, for both spin chan- 
nels. In these cases the system acts as an insulator in a particular energy region, i.e., it can act 
as a mutual spin switch. We have got conductance spectra as a function of energy of injected 
electrons both in weak- and strong-coupling limits. The sharp resonant peaks appear in the 
conductance spectrum in the weak-coupling limit, whereas in the strong-coupling limit, due to 
the broadening of energy levels, except for the presence of global gap, conductance does not 
show zero value. We have also presented the variation of conductance with the angle made by 
the local magnetic moment with the quantization direction, both including and excluding spin 
flip effect. 

Then, we have noticed the steplike behavior of I-V characteristics in the weak-coupling limit, 
which is a manifestation of quantized conductance of these kind of nanoscale systems. On the 
other hand, in the strong-coupling limit, the I-V spectra show almost continuous variation, 
instead of the step-like behavior. 

Throughout this presentation, we have made several important assumptions. Here we discuss 
briefly about these approximations. 

• All the calculations have been performed at OK, but the results will not significantly change 
in non-zero finite temperature region also, as long as the thermal energy (/c^T) is less than the 
average spacing energy levels of the magnetic conductor. 
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• Secondly we have neglected the electron-electron correlation effect. The inclusion of electron- 
electron correlation is a major challenge to us, since over the last few years a lot of efforts have 
been made to incorporate this effect, but no proper theory has yet been developed. 

• Electron-phonon interaction has also been neglected in our theoretical formulation, as we 
are dealing with OK temperature. But even at non-zero finite temperature also we could neglect 
this effect. Because the broadening of energy levels due to electron-phonon interaction is much 
smaller than that of due to coupling of the MQD device with source and drain. 

• The potential drops are considered at the boundaries only, not throughout the device. 

• All the results have been presented for ordered systems, but in presence of impurity the 
transport properties can change that we have not been studied yet. 

All the predicted results given above may be applied in designing spin polarized transistors 
in future. 
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